The moving barrier and an improved model for the multiplicity 

function 



A. Del Popolo 1 

Bogazigi University, Physics Department, 80815 Bebek, Istanbul, Turkey 

ABSTRACT 

I compare the numerical multiplicity function given in Yahagi, Nagashima & 
Yoshii (2004) with the theoretical multiplicity function obtained by means of the 
excursion set model and an improved version of the barrier shape obtained in 
Del Popolo & Gambera (1998), which implicitly takes account of total angular 
momentum acquired by the proto-structure during evolution and of a non-zero 
cosmological constant. I show that the multiplicity function obtained in the 
present paper, is in better agreement with Yahagi, Nagashima & Yoshii (2004) 
simulations than other previous models (Sheth & Tormen 1999; Sheth, Mo &: 
Tormen 2001; Sheth & Tormen 2002; Jenkins et al. 2001) and that differently 
from some previous multiplicity function models (Jenkins et al. 2001; Yahagi, 
Nagashima Sz Yoshii 2004) it was obtained from a sound theoretical background. 

Subject headings: cosmology: theory - large scale structure of Universe - galaxies: 
formation 



1. Introduction 



Two different kinds of methods are widely used for the study of the structure forma- 
tion. The first one is N-body simulations, that are able to follow the evolution of a large 
number of particles under the influence of the mutual gravity, from initial conditions to the 
present epoch. The second one are semi-analytical methods. Among them, Press-Schechter 
(hereafter PS) approach and its extensions (EPS) are of great interest since they allow us 
to compute mass functions (Press & Schechter 1974; Bond et al. 1991), to approximate 
merging histories (Lacey & Cole 1993, LC93 hereafter, Bower 1991, Sheth & Lemson 1999b) 
and to estimate the spatial clustering of dark matter haloes (Mo & White 1996; Catelan et 
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al. 1998, Sheth & Lemson 1999a). 

Although the analytical framework of the PS model has been greatly refined and ex- 
tended (Bond et al. 1991; Lacey and Cole 1993), it is well known that the PS mass function, 
while qualitatively correct, disagrees with the results of N-body simulations. In particular, 
the PS formula overestimates the abundance of haloes near the characteristic mass M* and 
underestimates the abundance in the high mass tail (Efstathiou et al. 1988; Lacey & Cole 
1994; Tozzi & Governato 1998; Gross et al. 1998; Governato et al. 1999). 

A better agreement between the numerical mass function and the analytic mass function 
can be obtained by incorporating into the PS ansatz the non-sphericity of collapse model 
(Del Popolo & Gambera 1998; Sheth & Tormen 1999 (hereafter ST) ; Sheth, Mo & Tormen 
2001 (hereafter SMT); Sheth & Tormen 2002 (hereafter ST1); Jenkins et al. 2001 (hereafter 
J01)), instead of the spherical model or taking into account the spatial correlation of density 
fluctuations (Nagashima 2001). 

More recently in order to investigate the functional form of the universal multiplicity 
function, Yahagi, Nagashima & Yoshii (2004) (hereafter YNY) performed five runs of N- 
body simulations with high mass resolution and compared them with different multiplicity 
function and with a fit by them proposed. 

They showed that discrepancies are observed between some of the quoted analytical 
multiplicity function with simulations. 

In the present paper, I shall use an improved version of the barrier shape obtained 
in Del Popolo & Gambera (1998), obtained from the parameterization of the nonlinear 
collapse discussed in that paper, taking account of asphericity and tidal interaction between 
protohaloes and the effects of a non-zero cosmological constant, together with the results of 
ST, ST1 in order to study the "unconditional" multiplicity function. 

The reasons that motivates this study are several: 

As previously reported, multiplicity functions like ST and J01, fit only approximatively 
high resolution N-body simulations like those of YNY, while the functional form proposed 
in YNY, provides a better fit when compared with the ST functional form. Unfortunately 
the functional form for the multiplicity function proposed in YNY and similarly that of J01 
(which is a fit to their "Hubble Volume" simulations of rCDM and ACDM cosmologies) are 
not based on any theoretical background. So it is important to find a better analytical form, 
which starting by "first principles" is able to fit in a better way simulations and is physically 
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motivated. I show that the function obtained in the present paper, similarly to that in YNY 
provides a better fit than the ST or other functional forms used in literature and moreover 
it has been obtained from solid physical, theoretical, arguments. 

The paper is organized as follows: in Sect. 2, 1 calculate the "unconditional" multiplicity 
function. Sect. 3 and 4 are devoted to results and to conclusions, respectively. 



2. The barrier model and the multiplicity function 



According to hierarchical scenarios of structure formation, a region collapses at time 
t if its overdensity at that time exceeds some threshold. The linear extrapolation of this 
threshold up to the present time is called a barrier, B. A likely form of this barrier is (ST, 
ST1): 



B(a 2 ,z) = v^[l + /3 (S/aSX] = V^cU 



1 + 



P 



[au 



(1) 



In the above equation a, (3 and a are constants, S* = 5 2 , where 5 c (t) is the linear extrapolation 
up to the present day of the initial overdensity of a spherically symmetric region, that 



collapsed at time t. Additionally, S = S* 



a{M) 



where a 2 (M) 



is the present day mass dispersion on comoving scale containing mass M. S depends on 
the assumed power spectrum. The spherical collapse model (SC) has a barrier that does 
not depend on the mass (eg. Lacey & Cole 1993 (LC93)). For this model the values of the 
parameters are a = 1 and (3 = 0. The ellipsoidal collapse model (EC) of ST has a barrier 
that depends on the mass (moving barrier). The values of the parameters are a = 0.707, 
f3 = 0.485, 7 = 0.615 and are adopted either from the dynamics of ellipsoidal collapse or 
from fits to the results of N-body simulations. 

In the following, I shall use an improved version of the barrier obtained in Del Popolo 
& Gambera (1998) to get the mass functions, which shall be compared with those obtained 
by PS, ST, J01, YNY, and with numerical simulations of YNY. Since the way the barrier is 
obtained is described in previous papers (see Del Popolo & Gambera 1998, 1999, 2000) the 
reader is referred to those papers for details. Assuming that the barrier is proportional to 
the threshold for the collapse, similarly to ST, the barrier can be expressed, in the case of a 
zero cosmo logical constant, in the form: 



B(M)=S c (u,z) = S c 



1 + 



r t J 2 ■ dr 
GMr 3 



(2) 



where <5 co = 1.68 is the critical threshold for a spherical model, r\ is the initial radius, r ta is 
the turn-around radius, I the specific angular momentum, a.\ = 0.585 and (3\ = 0.46. The 
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specific angular momentum appearing in Eq. (2) is the specific total angular momentum 
acquired by the proto-structure during evolution. In order to calculate L, I shall use the 
same model as described in Del Popolo & Gambera (1998, 1999) (more hints on the model 
and some of the model limits can be found in Del Popolo, Ercan & Gambera 2001, Sec. 3). 

Assuming a non-zero cosmological constant Eq. (2) is changed as follows (see Appedix): 



B(M) = 5 c (v,z) = 5 c 



rta r t J 2 ■ dr r ta r 2 



1+ / * o +A 



o GMr 3 6GM 



~ ft 
— f r 
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(3) 



where a 2 = 0.4 and f3 2 = 0.02 and Q A is the contribution to the density parameter coming 
from the cosmological constant. The values of ai, a 2 , f3\ and f3 2 are calculated so that the fit 
function at extreme right hand side of Eq. (3) reproduces the barrier shape (central part of 
Eq. (3) depending on I and A). ST1 connected the form of the barrier with the form of the 
multiplicity function. As shown by ST1, for a given barrier shape, B(S), the first crossing 
distribution is well approximated by: 

f(S)dS=\T(S)\eM-^)^ (4) 
where T(S) is the sum of the first few terms in the Taylor expansion of B(S): 

^£^§P (5) 

n=0 

The quantity Sf(S, t) is a function of the variable v alone. Since 8 C and a evolve with time in 
the same way, the quantity Sf(S,t) is independent on time. Setting 2Sf(S,t) = vf(v), one 
obtains the so-called multiplicity function f(u). The multiplicity function is the distribution 
of first crossings of a barrier B{y) by independent uncorrected Brownian random walks 
(Bond et al. 1991). That's why the shape of the barrier influences the form of the multiplicity 
function. 

In the excursion set approach, the average comoving number density of haloes of mass 
M the universal or "unconditional" mass function, n(M,z), is given by: 

(Bond et al. 1991), where p is the background density, In the case of the ellipsoidal barrier 
shape given in ST, namely Eq. 1 of the present paper, the Eqs. (4), (5), give, after truncating 
the expansion at n = 5 (see ST): 

vf(v) = vW 2 ^ 1 + (3(au 2 )- a g(a)} exp (-0.5az/ 2 [l + (3(au 2 )- a } 2 ) (7) 
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where 



g(a) =| 1 — a + 



a(a — 1) a(a — 1) • • • (a — 4) 

2! "' 5! 



(8) 



If the barrier takes account of the cosmological constant, like in Eq. (3), using the same 
method that lead to Eq. (7), we have that: 



uf(v) = A 1 1 + ai + a2 W — exp{- 
\ (az/) (az/) / V 27T 



■az/ 



. + 4 + a 



(az/) 1 (az/) 



Q2 



/2} (9) 



Using the values for (3 and a of ST (a = 0.707, 5 c (z) = 1.686(1 + ^), (3 ~ 0.485 and a ~ 0.615) 
in Eq. (7) we get (ST1): 



z//(z/) ~ A 2 1 + 



0.094 \ fcii/ 



1 + 



0.5 



az/ 



v0.6 



/2} 



(10) 



with A 2 ~ 1. This last result is in good agreement with the fit of the simulated first crossing 
distribution (ST): 

/ 1 \ I nu 

(11) 



vf(is)dis = A 3 ^1 + — yj ^ — exp(-az//2) 

where p = 0.3, and a = 0.707. 

The normalization factor A 3 has to satisfy the constraint: 



Jo 



fiy)dv = l 



(12) 



and as a consequence it is not an independent parameter, but is expressed in the form: 

A = [1 + 2- p 7r- 1/2 r(l/2 - p)] ^ = 0.3222 2 . (13) 

In the case of the barrier given in Eq. (2), the "unconditional" multiplicity function can 
be approximated by: 



z//(z/) ~ At 1 + 



(az/) 



0.585 




exp {—acu 



1 + 



d 



az/ 



v0.585 



} 



(14) 



where a = 0.707, b = 0.1218, c = 0.4019, d = 0.5526 and A 4 ~ 1.75 is obtained from the 
normalization condition. 



2 Note, that Eq. 11 gives a better fit to Eq. 7 if A ~ 0.3 and a ~ 0.79. Vice versa a smaller value of a 
(a ~ 0.63) and ^ = 1.08 in Eq. 7 gives a better fit to Eq. 11 (with A x = 0.3222 and a = 0.707), which was 
the one ST used to compare model and data. 
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In the case of the barrier with non-zero cosmological constant, Eq. (3), a good approx- 
imation to the multiplicity function is given by: 

2 



0.1218 0.0079\ lav 

} \0.585 + 7 \0^4 ) \/ 9T 

(au) (ais) J V Z7r 



vf(y) ~ A 5 [ 1 + ; ' g5 + —- ^ ) A / — exp {-0.4019az/ 



0.5526 0.02 

7 \0.585 + 7 70^4 



} 



(15) 

where A 5 — 1.75. As previously reported, for matter of completeness, to the previous 
functions, namely PS, ST, Eq. (15) we have to add J01, which satisfies the equation: 

uf(y) = 0.315exp(- | 0.61 + ln[a-\M)\ | 3 ' 8 ) (16) 

In order to express the above relation as a function of v, I substitute a~ l (M) = u/S c and 
I assume a constant value of <5 C , that of the Einstein-de Sitter Universe namely 5 C = 1.686. 
The above formula is valid for 0.5 < v < 4.8. 

YNY (Eq. 7, hereafter YNY7) proposed the following function to fit the numerical 
multiplicity function: 

uf(u) = A[l + (Bu/V2) c ]u D exp[-(Bu/V2) 2 }, (17) 

where, A is a normalization factor to satisfy the unity constraint, / °° f{y)dv = 1, therefore 

A = 2(B/V2) D {T[D/2] + F[(C + D)/2]}~ 1 . (18) 

The best-fit parameters are given as 5=0.893, C=1.39, and D=0.408, and from these pa- 
rameters, A is constrained so that A = 0.298. 

The CDM spectrum used in the present paper is that of Bardeen et al. (1986) (equa- 
tion (G3)). 



3. Results 

In this section, I compare the analytic multiplicity functions of PS, ST, J01, YNY7, and 
Eq. (15), of the present paper, with the numerical simulations of YNY. Those simulations 
adopt the ACDM cosmological parameters of Q m = 0.3, = 0.7, h = 0.7, and as = 1.0, 
using 512 3 particles in common (see YNY for details). 

The comparison between numerical multiplicity functions and theoretical ones is shown 
in Fig. 1. 

In the plot the solid line represents the multiplicity function obtained in the present 
paper, the short-dashed line YNY7, the dotted line the ST multiplicity function, the long- 
dashed line the J01 multiplicity function. The errorbars with open circles represents the run 
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140 of YNY, those with filled squares the case 70b, those with open squares the case 70a, 
those with filled circles the case 35b, those with crosses the case 35a. 

Since the data are available only in the region at v < 3, these functions could be 
erroneous at v > 3. 

Note that the comparison of the above curves, except for the PS model, with the results 
of N-body simulations show a very good agreement. However, there are some discrepancies 
between the YNY multiplicity function and other model functions (except this in the present 
paper). First, the multiplicity function of the present paper, similarly to that of YNY, in 
the \ow-u region of v < 1, systematically falls below the ST and the J01 functions. In this 
region the multiplicity function of the present paper is very close to that of YNY. 

As seen in Fig. 1, and in agreement with YNY, the numerical multiplicity functions 
reside between the ST and J01 multiplicity functions at 2 < v < 3 (except for the run 35b). 
Additionally, the numerical multiplicity functions have an apparent peak at v ~ 1 instead 
of the plateau that is seen in the J01 function. 

On the other hand, in the high-z/ region, where v is significantly larger than unity, 
the multiplicity function of the present paper like YNY takes values between ST and J01 
functions. These differences between numerical multiplicity functions and analytic ones, like 
ST, ST1 and J01, are within 1 a error bars, and they are possibly due to the different box 
sizes adopted (see YNY for a discussion). To be more precise, throughout the peak range of 
0.3 < v < 3, the ST multiplicity function is in disagreement with the high mass resolution N- 
body simulations of YNY and that of the present paper. As shown by YNY the ST functional 
form provides a good fit to them only choosing parameter values of a = 0.664, p = 0.321, 
and A 3 = 0.301. The multiplicity function obtained in the present paper has a peak at v ~ 1 
as in the ST function, and YNY numerical multiplicity function and YNY7, instead of a 
plateau as in the J01 function. 

I want to stress that the functional form proposed in YNY, namely YNY7, provides 
a better fit when compared with the ST functional form but it is not based on theoretical 
background. The function obtained in the present paper, similarly to YNY7 provides a better 
fit to simulations than the ST functional form, and at the same time has been obtained 
from solid physical, theoretical, arguments. The better agreement observed between the 
multiplicity function of the present paper and YNY simulations, when compared with the ST, 
is connected to the shape of the barrier (5 C ) . Taking account of the effects of asphericity and 
tidal interaction with neighbors, Del Popolo & Gambera (1998), showed that the threshold 
is mass dependent, and in particular that of the set of objects that collapse at the same time, 
the less massive ones must initially have been denser than the more massive, since the less 



- 8- 



massive ones would have had to hold themselves together against stronger tidal forces. 

The shape of the barrier given in Eq. (2) is a direct consequence of the angular momen- 
tum acquired by the proto-structure during evolution while Eq. (3) introduces the effects of 
the cosmological constant. 

Similarly to ST, the barrier increases with S (decrease with mass, M) differently from 
other models (see Monaco 1997a, b). It is interesting to note that the increase of the barrier 
with S has several important consequences and these models have a richer structure than 
the constant barrier model. 

The decrease of the barrier with mass means that, in order to form structure, more 
massive peaks must cross a lower threshold, S c (u,z), with respect to under-dense ones. At 
the same time, since the probability to find high peaks is larger in more dense regions, 
this means that, statistically, in order to form structure, peaks in more dense regions may 
have a lower value of the threshold, S c (v,z), with respect to those of under-dense regions. 
This is due to the fact that less massive objects are more influenced by external tides, 
and consequently they must be more overdense to collapse by a given time. In fact, the 
angular momentum acquired by a shell centred on a peak in the CDM density distribution is 
anti-correlated with density: high-density peaks acquire less angular momentum than low- 
density peaks (Hoffman 1986; Ryden 1988). A larger amount of angular momentum acquired 
by low-density peaks (with respect to the high-density ones) implies that these peaks can 
more easily resist gravitational collapse and consequently it is more difficult for them to 
form structure. Therefore, on small scales, where the shear is statistically greater, structures 
need, on average, a higher density contrast to collapse. 

It is evident that the effect of a non-zero cosmological constant adds to that of L. The 
effect of a non-zero cosmological constant is that of slightly changing the evolution of the 
multiplicity function with respect to open models with the same value of Qq. This is caused 
by the fact that in a flat universe with Q\ > 0, the density of the universe remains close 
to the critical value later in time, promoting perturbation growth at lower redshift. The 
evolution is more rapid for larger values (in absolute value) of the spectral index, n. 

As previously reported, the ST model gives a better fit to simulations than PS model, 
but it has some discrepancies with simulations. ST model was introduced at the beginning 
(Sheth & Tormen 1999) as a fit to the GIF simulations and in a subsequent paper (SMT) 
was recognized the importance of aspherical collapse in the functional form of the mass 
function. The effects of asphericity were taken into account by changing the functional 
form of the critical overdensity (barrier) by means of a simple intuitive parameterization 
of elliptical collapse of isolated spheroids. The model proposed in the present paper has 



- 9- 




Fig. 1. — The best-fit multiplicity function. In the plot the solid line represents the mul- 
tiplicity function obtained in the present paper, the short-dashed line YNY7, the dotted 
line the ST multiplicity function, the long-dashed line the J01 multiplicity function. The 
errorbars with open circles represents the run 140 of YNY, those with filled squares the case 
70b, those with open squares the case 70a, those with filled circles the case 35b, those with 
crosses the case 35a 
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Fig. 2. — Time dependence of the multiplicity function from the 35a run, for four redshift 
ranges of < z < 1 (open circles), 1 < z < 3 (open squares), 3 < z < 6 (open triangles), 
and z > 6, (crosses). Also shown are YNY7 (solid line) and the model of the present paper 
(dot-dashed line). 
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several similitudes with ST and ST1 models, namely it uses the excursion set approach as 
extended by ST1 to calculate the multiplicity function, but at the same time it differs from 
ST and ST1 for the way the barrier was calculated and for the fact that takes account 
of angular momentum acquisition, and a non-zero cosmological constant, things which are 
not taken into account into ST and ST1. These differences gives rise to a multiplicity 
function in better agreement with simulations. This shows the importance of the form 
of the barrier. The improvement of the model of the present paper and ST model with 
respect to PS is probably connected also to the fact that incorporating the non-spherical 
collapse with increasing barrier in the excursion set approach results in a model in which 
fragmentation and mergers may occur, effects important in structure formation. In the case 
of non-spherical collapse with increasing barrier, a small fraction of the mass in the Universe 
remains unbound, while for the spherical dynamics, at the given time, all the mass is bound 
up in collapsed objects. Moreover, incorporating the non-spherical collapse with increasing 
barrier in the excursion set approach results in a model in which fragmentation and mergers 
may occur (ST). If the barrier decreases with S (Monaco 1997a,b), this implies that all walks 
are guaranteed to cross it and so there is no fragmentation associated with this barrier shape. 

In other words, the excursion set approach with a barrier taking account effects of 
physics of structure formation gives rise to good approximations to the numerical multi- 
plicity function: the approximation goodness increases with a more improved form of the 
barrier (taking account more and more physical effects: angular momentum acquisition, 
non zero cosmological constant, etc). Another important aspect of the quoted method is 
its noteworthy versatility: for example it is very easy to take account of the presence of a 
non zero cosmological constant englobing it in the barrier. I recall that the YNY numerical 
multiplicity function assumes a non zero cosmological constant while the theoretical models 
(ST,ST1, J01) does not take this into account. 

Finally I checked the time dependence of the multiplicity function. Fig. 2 shows the 
multiplicity function from the 35a run, for four redshift ranges of < z < 1 (open circles), 
1 < z < 3 (open squares), 3 < z < 6 (open triangles), and z > 6, (crosses). At high redshifts, 
high-z/ halos in the exponential part of the YNY7 (solid line) function and Eq. (15) (dot- 
dashed line of the present paper) are probed. As redshift decreases, the probe window moves 
to the lower-z/ region. Fig. 2 shows that the multiplicity function of this paper, Eq. (15), 
and YNY7 both gives a good fit to the numerical simulations. For small values of z/, Eq. 
(15) is a slightly better fit to data, and at large values of v the two functions decays in the 
same way. 
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4. Conclusions 

In the present paper, I compared the numerical multiplicity function given in YNY with 
the theoretical multiplicity function obtained by means of the excursion set model and an 
improved version of the barrier shape obtained in Del Popolo & Gambera (1998), which 
implicitly takes account of tidal interactions between clusters and a non-zero cosmological 
constant. I showed that the barrier obtained in Del Popolo & Gambera (1998) gives rise 
to a better description of the multiplicity functions than other models (ST, J01) and the 
agreement is based on sound theoretical models and not on fitting to simulations. 

The main results of the paper can be summarized as follows: 

1) the non-constant barrier of the present paper combined with the ST1 model gives "un- 
conditional" multiplicity functions in better agreement with the N-body simulations of YNY 
than other previous models (ST, ST1, J01). 

2) The comparison of the theoretical multiplicity function of the present paper, in agreement 
with the YNY result, shows some discrepancies with the theoretical multiplicity function of 
several authors (ST, ST1, J01): e.g., the maximum value of the multiplicity function from 
simulations at v ~ 1 is smaller, and its low mass tail is shallower when compared with the 
ST multiplicity function. 

3) The multiplicity function of the present paper gives a good fit to simulations results as 
the fit function proposed by YNY, but differently from that it was obtained from a sound 
theoretical background. 

4) The excursion set model with a moving barrier is very versatile since it is very easy to 
introduce easily several physical effects in the calculation of the multiplicity function, just 
modifying the barrier. 

The above considerations show that it is possible to get accurate predictions for a number 
of statistical quantities associated with the formation and clustering of dark matter haloes 
by incorporating a non-spherical collapse which takes account of a non-zero cosmological 
constant in the excursion set approach. The improvement is probably connected also to the 
fact that incorporating the non-spherical collapse with increasing barrier in the excursion 
set approach results in a model in which fragmentation and mergers may occur, effects 
important in structure formation. Moreover, the effect of a non-zero cosmological constant 
adds to that of angular momentum slightly changing the evolution of the multiplicity function 
with respect to open models with the same value of matter density parameter. 
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5. Appendix 

The equation governing the collapse of a density perturbation taking account angular 
momentum acquisition by protostructures can be obtained using a model due to Peebles 
(Peebles 1993) (see also Del Popolo & Gambera 1998, 1999). 

Let's consider an ensemble of gravitationally growing mass concentrations and suppose that 
the material in each system collects within the same potential well with inward pointing 
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acceleration given by g(r) (see Del Popolo & Gambera 1998). We indicate with dP = 
f(L,rv r ,t)dLdv r dr the probability that a particle, of mass m, can be found in the proper 
radius range r, r + dr, in the radial velocity range v r — r, v r + dv r and with angular 
momentum L = mrvg in the range dL, or specific angular momentum I = L/m = rvg. The 
radial acceleration of the particle is: 



dv r 
~~dt 



- 9{r) 



l 2 (r) 



GM 



(19) 



where M is the mass of the central concentration. Eq. (19) can be derived from a potential 
and then from Liouville's theorem it follows that the distribution function, /, satisfies the 
collisionless Boltzmann equation: 



dt r dr dv r 



I 2 







(20) 



Assuming a non-zero cosmological constant Eq. (19) becomes: 

dv r GM l 2 (r) A 
— - = 1 — H r 

dt r 2 r 3 3 



(21) 



(Peebles 1993; Bartlett & Silk 1993; Lahav 1991; Del Popolo & Gambera 1998, 1999). 
Integrating Eq. (21) we have: 



1 f dr 
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12 j A 2 
— dr + — r + e 



(22) 



where the value of the specific binding energy of the shell, e, can be obtained using the 
condition for turn- around, ^ = 0. 

In turn the binding energy of a growing mode solution is uniquely given by the linear 
overdensity, 5 iy at time ij. From this overdensity, using the linear theory, we may obtain 
that of the turn-around epoch and then that of the collapse. We find the binding energy of 
the shell, C, using the relation between v and Si for the growing mode (Peebles 1980) in Eq. 
(22) and finally the linear overdensity at the time of collapse is given by: 



1 + 



rta n j2 . dr Tjy_ 

GMr 3 QGM 



(23) 



where a x = 0.585, fa = 0.46, a 2 = 0.4 and (3 2 = 0.02 



